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Abstract 

This paper focusses attention on the strange nonchaotic attractors (SNA) 
of a quasiperiodically forced dynamical system. Several routes, including the 
standard ones by which the appearance of strange nonchaotic attractors takes 
place, are shown to be realizable in the same model over a two parameters 
(/-e) domain of the system. In particular, the transition through torus dou- 
\Q • bling to chaos via SNA, torus breaking to chaos via SNA and period doubling 

bifurcations of fractal torus are demonstrated with the aid of the two pa- 

00 , 

\ rameter (f-e) phase diagram. More interestingly, in order to approach the 

strange nonchaotic attractor, the existence of several new bifurcations on the 
torus corresponding to the novel phenomenon of torus bubbling are described. 
Particularly, we point out the new routes to chaos, namely, (1) two frequency 
quasiperiodicity — > torus doubling — > torus merging followed by the gradual 
c~| | fractalization of torus to chaos, (2) two frequency quasiperiodicity — > torus 

J-^ ' doubling —* wrinkling — > SNA — > chaos — > SNA — > wrinkling — > inverse torus 

doubling — > torus — > torus bubbles followed by the onset of torus breaking to 
^ ■ chaos via SNA or followed by the onset of torus doubling route to chaos via 

SNA. The existence of the strange nonchaotic attractor is confirmed by calcu- 
lating several characterizing quantities such as Lyapunov exponents, winding 
numbers, power spectral measures and dimensions. The mechanism behind 
the various bifurcations are also briefly discussed. 

PACS number(s): 05.45. +b 
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I. INTRODUCTION 



In nonlinear dynamical systems strange nonchaotic attractors (SNA) are considered as 
complicated structures in phase space, which is a property usually associated with chaotic 
attractors. The pioneering work of Grebogi et al. [1] revealed that there are some possi- 
bilities of strange attractors in certain types of dynamical systems which are not chaotic. 
These strange attractors are strange in the spirit that geometrically they are strange (fractal 
dimensional) objects in phase space. On the other hand, they would not exhibit sensitivity 
to initial conditions (for example, Lyapunov exponents are negative), and hence are not 
chaotic. These strange nonchaotic attractors can arise in physically relevant situations such 
as quasiperiodically forced pendulum [2-4], quantum particles in quasiperiodic potentials 
[5], biological oscillators [6], Duffing type oscillators [7-10], velocity dependent oscillators 
[11], electronic circuits [12,13] and in certain maps [14-22]. Also these exotic attractors were 
confirmed by an experiment consisting of a quasiperiodically forced, buckled magnetoelastic 
ribbon [23], in analog simulations of a multistable potential [24] and in a neon glow discharge 
experiment [25]. 

While the existence of strange nonchaotic attractors has been firmly established, a ques- 
tion that remains interesting is what are the possible routes by which they arise and ul- 
timately become chaotic, and how do these attractors are born in a system (mechanism). 
Several routes have been identified in recent times and for a few of them typical mechanisms 
have also been found for the creation of SNA. The major routes by which the SNA appear, 
may be broadly classified as follows. 

1. Ding et al. route [6]: two frequency quasiperiodicity — > three frequency quasiperi- 
odicity — > strange nonchaotic attractor — > chaos. 

2. Kapitaniak et al. route [10]: two frequency quasiperiodicity — > strange nonchaotic 
attractor — > three frequency quasiperiodicity — > chaos. 

3. Heagy and Hammel route [14]:two frequency quasiperiodicity — > torus doubling 
— > wrinkling — > strange nonchaotic attractors — > chaos. 

4. Feudal et al. route [17]: two frequency quasiperiodicity -^wrinkling strange 
nonchaotic attractors — > chaos. 

5. Yalencinkaya and Lai route [9]: two frequency quasiperiodicity — > strange non- 
chaotic attractor (on-off intermittency type attractor) — > chaos. 

6. Venkatesan and Lakshmanan route [11]: two frequency quasiperiodicity — > 
torus doubling — > torus merging — > wrinkling — > strange nonchaotic attractor — > chaos. 

7. Kapitaniak and Chua route [13]: two frequency quasiperiodicity — > strange 
nonchaotic trajectories on torus — > chaos. 

8. Nishikawa and Kaneko route [21]:two frequency quasiperiodicity — > strange 
nonchaotic attractors — > chaos. 

Different mechanisms have been identified for some of the above routes. In particular, it 
has been shown that the birth of SNA in the Heagy and Hammel [14] route is due to the 
collision between a period doubled torus and its unstable torus. Feudal et al. [17] explained 
in their route that the SNA also appears a result of a collision of stable and unstable torus 
in a dense of set of points. However, Nishikawa and Kaneko [21] discussed in their route 
that the SNA emerges without interaction of stable and unstable torus. Moreover, the loss 
of transverse stability of the torus can also lead to the birth of SNA, as in the case of 
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the Yalencinkaya and Lai route [9] above. For other routes as far as our knowledge goes 
mechanisms have not yet been found. 

Also, most studies of strange nonchaotic attractors have focussed on their characteri- 
zation using the spectral properties [2, 5], geometrical properties [6], local divergence of tra- 
jectories[10], phase sensitivity and rational bifurcations [15-18], and functional maps and 
invariant curves [20-21]. 

In this paper, we demonstrate the existence of at least five different routes to chaos via 
strange nonchaotic attractors in a single dynamical system namely a quasiperiodically forced 
velocity dependent nonpolynomial oscillator system over a two parameters (/-e) space. To 
start with the birth of the strange nonchaotic attractors associated with two important 
routes, namely, (i) torus breaking and (ii) torus doubling, have been studied in our model. 
In low dimensions, Bier and Bountis have shown that a dynamical system that undergoes 
one or more period doublings need not complete the entire infinite Feigenbaum cascade, but 
it may be possible to have only a finite number of period doublings, followed by for example 
undoublings or other bifurcations [28]. The possibility of such a novel remerging bifucation 
phenomenon of the torus doubling sequence in the quasiperiodically forced system has not yet 
been reported. Since the system that we consider possesses more than one control parameter 
and remains invariant under the reflection symmetry, the remergence is likely to occur as 
in the case of low dimensional systems [28,29]. To confirm such a possibility, our numerical 
studies show that in some regions of the (f-e) parameter space, torus doubled orbit emerges 
and remerges from a single torus orbit at two different parameter values of e to form a torus 
bubble. Such a remerging bifurcation can lead to tame the growth of the torus doubled trees 
and the development of the associated universal route to chaos further. However, the nature 
of remerging torus doubled trees or more specifially torus bubbling, ensures the existence 
of different routes for the creation of SNA when the full range of parameters are taken into 
account. To illustrate these possiblities in our system, we enumerate two new types of routes 
as: (1) two frequency quasiperiodicity — > torus doubling — > torus merging followed by the 
gradual fractalization of torus to chaos [11], (2) two frequency quasiperiodicity — > torus 
doubling — > wrinkling SNA — > chaos — > SNA — > wrinkling — > inverse torus doubling — > 
torus — > torus bubbles followed by the onset of torus breaking to chaos via SNA or followed 
by the onset of torus doubling route to chaos via strange nonchaotic attractor. Finally, we 
also show the occurrence of period doubling bifurcations of the destroyed torus (strange 
nonchaotic attractor) in our model. 

This paper is organised as follows. Sec. 2 describes the system and the salient fea- 
tures of its dynamics. Sec. 3 describes some of the characterizing quantities of strange 
nonchaotic attractors in comparison with chaotic attractors such as Lyapunov exponents, 
winding number, power spectral analysis and dimensions. These quantities have been used 
to distinguish between quasiperiodic, strange nonchaotic and chaotic attractors. The birth 
of strange nonchaotic attractors from the transitions of the two frequency quasiperiodic at- 
tractors are observed in the five different routes mentioned above in sec. 4. The first one is 
that of torus breaking to chaos via SNA. The second one is through torus bubbling followed 
by the gradual fractalization to chaos. The third route is the transition from torus doubling 
to chaos via strange nonchaotic attractor. The next one we notice is the possibilty of torus 
doubling to chaos via SNA followed by the inversely advancing type of torus. The last one 
we observe is the period doubling bifurcation of the destroyed torus. Finally, in sec. 5 we 
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summarize our results. 



II. QUASIPERIODICALLY FORCED VELOCITY DEPENDENT SYSTEM 

Let us consider briefly the dynamics of a damped and driven rotating parabola system 
and discuss some of its general properties as reported in ref. [11]. 

A mechanical model describing the motion of a freely sliding particle of unit mass on a 
parabolic wire (z=\/Xx 2 ) rotating with a constant angular velocity f2 (f2 2 = Vl 2 = —u 2 + 
9\f{\) and g is the acceleration due to gravity)can be associated with a velocity dependent 
Lagrangian [27] 



L = - 

2 



1 + Xx 2 )x 2 - uj 2 Q x 2 . (1) 



Here overdot stands for a derivative with respect to time. The corresponding equation of 
motion is 

(1 + Xx 2 )x + Xxx 2 + u 2 x = 0, (2) 

When uj 2 > 0, Eq.(2) can be integrated in terms of elliptic integrals. Interesting bifurca- 
tions and different routes to chaos occur in the above model when the system is acted upon 
by additional damping and external forcing [11]. In this case, Eq.(2) gets modified to 

(1 + Xx 2 )x + Xxx 2 + uj 2 x + ax = f cos out. (3) 

The familiar period doubling bifurcations, preceded by a symmetry breaking bifurcation, 
intermittency and antimonotonicity have been identified by us earlier in ref [11]. 

Another interesting physical situation is the case in which there is an additional para- 
metric modulation in the angular velocity, 

Q = o (l + e cos ujpt), 

so that we can replace oo 2 = gyfX — Q 2 in Eqs. (l)-(3) by gy/\ — f22 = uo 2 — Vl 2 \2ecosuJ p t + 
0.5e 2 (l + cos2uj p t)]. Then the equation of motion becomes (see for example p. 351 in ref. [27] 
and also ref. [11]) 

(1 + Xx 2 )x + Xxx 2 + uj 2 x — £l 2 D [2e cos u p t + 0.5e 2 (l + cos2uj p t)]x + ax = f cosuj e t, (4) 

where e is a small parameter. 

We have already noted an interesting new quasiperiodic route to chaos, namely two 
frequency quasiperiodicity — > torus doubling — > torus merging — > wrinkling — > strange non- 
chaotic attractor — > chaos in the system (4) in ref. [11]. In this paper we make a detailed 
study of this and many other quasiperiodic routes to chaos which can occur in this model in 
a range of f-e parameter values and compare them. For our analysis we rewrite the system 
(4) as 
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X 



-Xxy 2 - {uj 2 - tt 2 Q [2e cos <p + 0.5e 2 ( 1 + cos 20)] }x - ay + f cos 



[l + Xx 2 ) 



4> = Up, 
9 = u e . 



(5) 



We note that the system (5) remains invariant under the reflection symmetry (x, y, f) 
— > (—x,—y,—f), (or equivalently (4) under the transformation (x,f) — > (— x, — /)). In 
analogy with low dimensional systems involving more than one control parameter when 
period bubbles occur [28], one may expect remergence of torus doubling sequences to occur 
in this model, which we indeed show to be true in the following. 



III. CHARACTERIZATION OF THE QUASIPERIODIC, STRANGE 
NONCHAOTIC AND CHAOTIC ATTRACTORS 

There are several quantities to characterize the attractors, which are useful to distinguish 
strange nonchaotic from chaotic and quasiperiodic attractors. We briefly review some of 
them which we will use in our study. 



A. The Lyapunov exponents 

For the system (5), there are two Lyapunov exponents that are trivial in the sense that 
they are identically zero by virtue of the two excitation frequencies. Let the Lyapunov 
exponents \ be ordered by their values, Ai > A2 > A3 > A4. We then have the following 
possibilities: (1) two frequency quasiperiodic attractors, Ai = A2 = > A3,Azi, (2) three 
frequency quasiperiodic attractors, Ai = A2 = A3 = > A4, (3) strange chaotic attractors, 
atleast Ai > and (4) strange nonchaotic attractors, the same as two frequency quasiperiodic 
attractors. 



B. Winding number 

The winding number for the orbit of Eq.(5) is defined by the limit 

iy=lim{ (a(t) - Q(t " )) }, (6) 

where (x,x) = (r cos a, r sin a). For the two frequency quasiperiodic attractors, the winding 
number satisfies 

W = -u p + —u e , 7 
n n 

where l,m, and n are integers. Combining the winding number and the Lyapunov exponents, 
we can distinguish the strange nonchaotic attractors from the other nonchaotic attractors 
as noted in Table I. 
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C. Power spectrum analysis 



To quantify the changes in the power spectrum (obtained using Fast Fourier Tranform 
(FFT) technique) one can compute the so called spectral distribution function N(a), de- 
fined to be the number of peaks in the Fourier amplitude spectrum larger than some 
value say a. Scaling relations have been predicted for N(a) in the case of two, three 
frequency quasiperiodic and strange nonchaotic attractors. These scaling relations are 
7V"(cr) ~ \n^,N(a) ~ In 2 a, and N(a) ~ a^ 13 respectively corresponding to two, three 
frequency quasiperiodic and strange nonchaotic attractors. In the Romeiras and Ott [2] 
studies, the power law exponent was found empirically to lie with in the range 1 < (5 < 2 for 
the strange nonchaotic attractor. Thus the above characterizations allow us to distinguish 
the strange nonchaotic attractor from other nonchaotic attractors as seen in Table I. 



D. Dimensions 

To quantify geomtric properties of attractors, several methods have been used to com- 
pute the dimension of the attractors. Among them, what we have used is the correlation 
dimension (introduced by Grassberger and Procaccia [26]) which may be computed from 
the correlation function C(R) defined as 



C(R) = lim 

N^oo 



N 



N 2 



H(R- \ Xi - Xj \) 



where Xi and Xj are points on the attractor, H(y) is the Heaviside function (1 if y > and 
if y < 0), N is the number of points randomly chosen from the entire data set. The Heaviside 
function simply counts the number of points within the radius R of the point denoted by Xi 
and C(R) gives the average fraction of points. Now the correlation dimension is defined by 
the variations of C(R) with R: 



C{R) ~ R d 



as 



R->0. 



Therefore the correlation dimension (d) is the slope of a graph of log C(R) versus log R. Once 
one obtains the dimensions of the attractors, it will be easy to quantify strange property of 
the attractors as seen in Table I. 



IV. DIFFERENT ROUTES TO CHAOS VIA STRANGE NONCHAOTIC 

ATTRACTORS 

Now we consider the combined effect of both the external and parametric forcings in 
Eq. (5). To be concrete, we consider the dynamics of (5) and numerically integrate it 
using the fourth order Runge-Kutta algorithm with adaptive step size with the values of 
the parameters fixed at u; 2 =0.25, A=0.5, a=0.2, i7 2 =6.7, lo p =1.0 and cj e =0.991. Various 
characteristic quantities such as the winding numbers, Lyapunov exponents, power spectral 
measures and dimensions as discussed in the previous section have been used to distinguish 
quasiperiodic, strange nonchaotic and chaotic attractors. Further to identify the different 
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attractors, the dynamical transitions are traced out by two scanning procedures: (i) varying 
/ at a fixed e and (ii) varying e at a fixed /. The resulting phase diagram in the (/-e) 
parameter space is shown in Fig.l. The various features indicated in the phase diagram are 
summarised and the dynamical transitions are elucidated in the following. 

A. Torus breaking bifurcations and the birth of strange nonchaotic attractors 

For low / and low e values, the system exhibits two frequency quasiperiodic oscillations 
denoted by IT in Fig. 1. When the value of e exceeds a certain critical value for a fixed low 
/, a transition from two frequency quasiperiodic (IT) to chaotic attractor(C) via strange 
nonchaotic attractor (S) occurs on increasing e. For example, we fix the strength of the 
external forcing parameter value as /=0.302 and vary the modulation parameter e. For 
e=0.03, Fig. 2a of the attractor has smooth branches and this indicates that the system 
is in a two frequency quasiperiodic state. As e increases, the branches in the Fig. 2b 
start to wrinkle (Wl). As e increases further, the attractor becomes extremely wrinkled 
and has several sharp bends. The sharp bends appear to become actual discontinuties at 
e=0.0419 and ultimately result in fractal phenomenon. Such a phenomenon is essentially 
the result of the collision of stable and unstable torus in a dense set of points as was shown 
by Feudal et al. in their route to chaos via SNA [17]. At such values, the nature of the 
attractor is strange (Fig. 2c) eventhough the largest Lyapunov exponent in Fig. 3 remains 
negative. For this attractor, the correlation dimension is 1.33 while the Fourier amplitude 
scaling constant (a) is 1.54. Winding number W does not satisfy the realtion (7) for this 
attractor. Hence, these studies confirm further that the attractor shown in Fig. 2c is strange 
nonchaotic. As e increases further, an attractor visibly similar to Fig. 2c appears (see Fig 
2d for /=0.042). However, it has a positive Lyapunov exponent and hence it corresponds 
to a chaotic attractor. 

B. Remerging torus doubling bifurcations: torus bubble and its consequences 

1. Torus bubbling 

On increasing the forcing parameter / further, 0.305< / <0.325, the fascinating novel 
phenomenon of torus bubble appears within a range of values of e. Within this range of /, 
on increasing the value of e along the same line, the onset of chaos is realised via strange 
nonchaotic attractor. To be more specific, the parameter / is fixed at 0.32 and e is varied. 
For e=0.03, the attractor is a two frequency quasiperiodic attractor (Figs. 4a & 5a). As e 
is increased to e=0.0317, the attractor undergoes a torus doubling bifurcation (Figs. 4b & 
5b). The corresponding period doubled torus attractor is denoted as 2T in Fig. 1. We note 
from Figs. 4 and 5 that the two strands in the (x, (p) projection become four strands when 
torus doubling bifurcation occurs. When we compute <p modulo 4n instead of 2n during 
integration, we notice from Figs. 5 that the two bifurcated strands of length 2ir are actually 
a single strand of length 4n. As a result, it can be concluded that the torus doubling is 
nonetheless a length doubling bifurcation. Further, it may noted that this bifurcation is 
geometrically very similar to that of period doubling bifurcation in the three dimensional 
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flows. One then expects as e is increased further that the doubled attractor has to continue 
the doubling sequence as in the case of period doubling phenomenon. Instead, in the present 
case, interestingly the strands of the length doubled attractor begin to merge into that of a 
single attractor at e=0.0353 as shown in Fig 6a, leading to the formation of a torus bubble 
(see Fig. 1), reminiscent of period bubbles in low dimensional systems [28,29]. On further 
increase of the value of e, the transition from two frequency quasiperiodicity to chaos via 
strange nonchaotic attractor takes place due to torus breaking bifurcations as discussed in 
section 4. A (see Figs 6 b,c,& d and Fig. 7). 

It has been argued in the case of period bubbling in low dimensional systems [28] that 
the cause of formation of the period bubbles is essentilally due to the presence of reflection 
symmetry combined with more than one control parameter present in the system. It appears 
that similar argument holds good for the case of higher dimensions for the formation of torus 
bubbles. 

2. Formation of multibubbles 

As the forcing parameter / is increased further in the region 0.325< / <0.332, the 
evolution of attractor undergoes the following transition to chaos, wherein more than one 
bubble is formed on increasing the value of e: two frequency quasiperiodicity — > torus 
doubling — > wrinkling — > inverse torus doubling (doubled torus) — > merged torus — > torus 
bubble — > merged torus — » wrinkling — > SNA — > chaos. To illustrate this possibility, let us 
fix the forcing parameter value as /=0.328 and vary the e value. For e=0.03 the attractor is 
two frequency quasiperiodic. As e is increased to e=0.0313, the attractor undergoes to torus 
doubling bifuraction. The doubled attractor begins to wrinkle when the e value is increased. 
However, this wrinkled attractor appears to become again a torus doubled attractor, instead 
of approaching the SNA while the e value is increased further. This doubled attractor merges 
into a single torus through inverse bifurcation on increasing the value of e. The merged torus 
again forms a torus bubble and then finally transits to chaos via wrinkling and SNA as the 
value of the e is incresed further. 

On incresing the forcing parameter / further, 0.332< / <0.335, the transition from 
two frequency quasiperiodicity to chaos via SNA takes place through the following route, 
wherein more than two bubbles are formed as e increases:two frequency quasiperiodicity — > 
torus doubling — > wrinkling — > inverse torus doubling — > merged torus — > torus bubble — ► 
torus — > torus bubble — > torus — > wrinkling — > SNA — > chaos. 

C. Strange nonchaotic and chaotic attractors within and outside the main torus 

bubble 

On further increase of /, / >0.335, inside the main torus bubble we observe interesting 
possiblities of strange nonchaotic and chaotic attractors via wrinkling as e increases. Then 
two interesting possibilities arise inside the main bubble. The dynamics outside the main 
bubble more or less follows the previous case 3B. The details are as follows. 
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1. SNA within the main torus bubble 



In a narrow region of /, 0.335< / <0.339, the SNA undergoes an inverse bifurcation 
scheme leading to two frequency quasiperiodic attractor as e increases through the following 
route: two frequency quasiperiodicity — > torus doubling — > wrinkling — >SNA — > wrinkling 
— > inverse torus doubling (doubled torus) — > merged torus. For example, the forcing pa- 
rameter / is fixed at /=0.337 and e is varied. For e=0.03 the attractor is a two frequency 
quasiperiodic one (Fig. 8a ). As e is increased to e=0.031, the attractor undergoes a torus 
doubling bifuraction ( as seen in Fig. 8b). In lower dimensional systems, the period doubling 
occurs in an infinite sequence until the accumulation point is reached, beyond which chaotic 
behaviour appears. However, with tori, in the present case, the truncation of the torus dou- 
bling begins when the two strands become extremely wrinkled (W2) when the e value is 
increased, as shown in Fig. 8c These strands lose their continuity as well as smoothness and 
become strange at e =0.0339. At such values, the attractor possesses geometrically strange 
property but does not obey the sensitivity to the initial conditions (the maximal Lyapunov 
exponent is negative as seen in Fig. 9) and so it is named as strange nonchaotic attractor 
(Fig. 8d). The emergence of such SNA is due to the collision of stable doubled torus and 
its unstable parent as was shown by Heagy and Hammel[14] in their route. Interestingly the 
SNA (Fig. 8e), instead of approaching a chaotic attrator as the e value increases, becomes 
wrinkled (Fig. 8f) and then torus doubled attractor (Fig. 8g). The doubled attractor again 
merges into a single torus (Fig. 8h) on further increasing the value of e. 

2. Chaotic attractor within the main torus bubble 

In a rather large region of /, / >0.339, the SNA as formed above transits into chaotic 
attractor on increasing the value of e further through the following route:two frequency 
quasiperiodicity — > torus doubling — > wrinkling — > SNA — > chaos — > SNA — > wrinkling 
— > inverse torus doubling (doubled torus) — > merged torus. To illustrate this possibility, 
let us choose the parameter /=0.342, and vary the value of e. For e=0.03 the attractor 
is a two frequency quasiperiodic attractor (Fig. 10a). As e is increased to e=0.0309, the 
attractor undergoes a torus doubling bifurcation ( Fig. 10b ). As e is increased further 
the strands of the doubled attractor begin to wrinkle (W2), as shown in Fig. 10c. The 
formation of sharp bends in the strand of the attractor is now clear as e is increased further. 
These bends tend to become actual discontinuities at e= 0.0337, as shown in Fig. lOd. 
The emergence of such dicontinuities on the torus is due to the collision of stable doubled 
torus and its unstable parent which is similar to the one found by Heagy and Hammel 
[14]. At such values, the attractor loses smoothness and becomes "strange". The attractor 
shown in Fig. lOd is nothing but strange nonchaotic as the maximum Lyapunov exponent 
works out to be A=— 0.01213 (Fig. 11). Further the correlation dimension is 1.49, the scaling 
constant (a) is 1.38 and winding number W does not satisfy the relation (7) for this attractor. 
Hence, these characteristic studies confirm further that the attractor shown in Fig. lOd is 
strange but nonchaotic. On increasing e value further, for a narrow range of e values, period 
doubling bifurcation of SNA is also noticed, which will be discussed in detail in the following 
subsection. On further increase of the value of e to 0.034, we find the emergence of a chaotic 
attractor (Fig. lOe), which though visibly similar to the nonchaotic strange attractor Fig. 
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lOd, has a positive Lyapunov exponent (see Fig. 11). The chaotic attractor again becomes 
SNA when e is further increased (Fig. lOf). As the value of e is still increased, the SNA 
becomes torus doubled attractor (Fig. lOg) via wrinkling. This doubled attractor then 
merges into a single torus (Fig. lOh) when the value e is continuously increased. 

3. Dynamics outside the main torus bubble 

Two further interesting transitions exist outside of the main torus bubble, namely (i) 
torus breaking to chaos via SNA and (ii) torus doubling to chaos via SNA. The details are 
as follows. 

In a narrow region of /, 0.335< / <0.345, the transition from two frequency quasiperi- 
odicity to chaos via SNA takes place outside the main bubble through the following route 
as e increases beyond e =0.0361: torus — > torus bubble — * torus — > torus bubble — > torus — > 
wrinkling — > SNA — > chaos. However in the region of /, 0.345< / <0.352, one also observes 
a transition from a wrinkled two torus (W2) to a wrinkled one torus (Wl). 

Further increase of the value of / beyond 0.352, for e values greater than 0.0362, intro- 
duces yet another kind of transitions beyond the main bubble as torus — > torus doubling 
— > wrinkling — > inverse torus doubling — > merged torus — > torus doubling — > wrinkling — > 
strange nonchaotic attractor — > chaos. 

D. Period doubling bifurcations of destroyed torus (SNA) within the main torus 

bubble 

In the previous subsections, we have seen that the period doubling bifurcation of torus 
has been truncated by the destruction of the torus leading to the emergence of strange non- 
chaotic attractor in certain regions of the (f-e) parameter space. However, we observe in the 
present system that in some cross sections of the (f-e) parameter space, the period doubling 
bifurcation phenomena still persists in the destroyed torus, eventhough the actual doubling 
sequence of the torus has been terminated. Such a route has also been recently observed in 
coupled Duffing oscillators [8] and in certain maps [22]. The doubling of destroyed torus has 
been observed in the present model in a rather long range of /, 0.338< / <0.358 and for a 
narrow range of e values denoted by S2 in Fig.l. For example, let us choose /=0.345 and 
vary the value of e. For e=0.03, the attractor is a two frequency quasiperiodic torus (Fig. 
12a). As e is increased to e=0.0305,the system undergoes torus doubling bifurcations (Fig. 
12b). On increase of the value of e further, the attractor begins to wrinkle and finally ends 
up with fractal nature (SNA) (Fig. 12c). On further increase of the value of e, the fractal 
torus undergoes doubling bifurcation (Fig. 12d). If the parameter e is still increased, the 
doubled fractal torus merges into a single fractal torus and finally transits into the chaotic 
attractor. 

V. CONCLUSION 

In this paper, we considered the dynamics of the specific example of the quasiperiodically 
forced velocity dependent nonpolynomial oscillator system (5) which illustrates many of the 
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typical routes to chaos via strange nonchaotic attractors. It was found that the first two of 
these routes can be realized in the following ways. 

• two frequency quasiperiodicity — » torus doubling — > wrinkling — > strange non- 
chaotic attractor — > chaos. Here, the birth of SNA is due to the collision of stable 
doubled orbit and its unstable parent torus [14] . 

• two frequency qusiperiodicity ^wrinkling — > strange nonchaotic attractor — > 
chaos. In this case, the emergence of SNA is essentially the result of interac- 
tion of stable and unstable torus in a dense set of points [17]. 

More interestingly, we have pointed out the novel possibility of the torus bubbling. That 
is torus doubling bifurcations in dynamical systems can under suitable circumstances form 
finite sequence which 'merge' in some cross sections of the parameters space, inhibiting the 
onset of torus doubling route to chaos. Such remerging bifurcations having finite number of 
'bubbles' occur only within some range of the parameters values. An important consequence 
is of such remerging is that the orbits become again stable and relatively large regions 
reappear around them, where the motion is regular and predictable. To illustrate such 
remerging torus doubling bifurcations, in our present study, we have shown two more routes: 

• two frequency quasiperiodicity — > torus doubling — > torus merging followed by 
the gradual fractalization of torus to chaos via SNA 

• two frequency quasiperiodicity — > torus doubling — > wrinkling — > SNA — > chaos 
— > SNA — > wrinkling — > inverse torus doubling — > torus — > torus bubbles followed 
by the onset of torus breaking to chaos via SNA or followed by the onset of torus 
doubling route to chaos via SNA. 

From these routes, it can be concluded that prior to standard routes for transition to strange 
nonchaotic attractor, the possibilities of several bifurcations on the torus can be realised. 

Finally, the period doubling bifurcations of the destroyed torus have also been observed 
in our model in a narrow region of the (/ — e) parameter space: 

• two frequency quasiperiodicity — > torus doubling — > wrinkling — > destroyed torus 
— > period doubling of destroyed torus — > merged destroyed torus — > chaos. 
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TABLES 

TABLE I. Characteristics of attractors 



Types of attractors Winding number Lyapunov exponents Power law relations Dimensions 

three frequency W / ^uj p + Ai < 0, A 2 = A 3 = A 4 = N(a) = In 2 a Integers 
quasiperiodic 

two frequency W =y 1 uj v + n We Ai, A2 < 0, A3 = A4 = N(a) = In ^ Integers 
quasiperiodic 

strange nonchaotic W / ^uj p + ^ui e Ai, A2 < 0, A3 = A4 = N(a) = a' 13 Fractals 

attractor 1 < (3 < 2 
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FIGURES 



FIG. 1. Phase diagram of the two parameters (/ — e) space exhibited by the system (5). 
Regions of different attractors are denoted as: IT - two frequency quasiperiodicity attractor, 2T- 
torus doubled attractor, Wl- wrinkled attractor of period one, W2 - wrinkled attractor of period 
two, S - strange nonchaotic attractor, S2- period doubled strange nonchaotic attractor and C - 
chaotic attractor. 

FIG. 2. Projection of the two frequency quasiperiodic attractors of Eq. (5) for /=0.302: 
Poincare plot with mod 2ir in the (x,^>) plane, (a) two frequency quasiperiodic attractor at 
e=0.030, (b) torus wrinkled attractor for e=0.0405, (c) strange nonchaotic attractor for e=0.0419, 
(d) chaotic attractor for e=0.042, The other parameters are 0^=0.25, A=0.5, a=0.2, u> p =1.0, 
n 2 Q =6.7 and w e =0.991. 

FIG. 3. Largest Lyapunov exponent X ma x vs.e corresponding to fig. 2 

FIG. 4. Projection of the two frequency quasiperiodic attractors of Eq. (5) for /=0.302: 
Poincare surface of section in the (x,y) plane (a) torus at e=0.03, (b)torus doubled attractor 
at e=0.0317. The other parameters are ^=0.25, A=0.5, a=0.2, uj p =1.0, f^=6.7 and u e =0.991. 

FIG. 5. Projection of the two frequency quasiperiodic attractors of Eq.(5) for 7=0.302: 
Poincare surface of section with (f> mod 2ir in the (x,</>) plane (a) two frequency quasiperiodic 
attractor for e=0.03, (b) torus doubled attractor at e=0.0317, (c) & (d) same as (a) & (b) except 
(f) mod 47r during integration. The other parameters are Wq=0.25, A=0.5, a=0.2, u p =1.0, 0,^=6.7 
and w e =0.991. 

FIG. 6. Projection of the two frequency quasiperiodic attractors of Eq.(5) for /=0.32: Poincare 
plot with <p mod 2tt in the (x,^>) plane, (a) merged attractor for e=0.0353 (b) wrinkled attractor 
for e=0.0399 (c) strange nonchaotic attractor for e=0.041 (d) chaotic attractor for e=0.0413 The 
other parameters are 0^=0.25, A=0.5, a=0.2, u) p =1.0, ^=6.7 and cj e =0.991. 

FIG. 7. Largest Lyapunov exponent \ ma x vs.e corresponding to Figs. 4,5 & 6 

FIG. 8. Projection of the two frequency quasiperiodic attractors of Eq.(5) for /=0.337: 
Poincare plot with (f> mod 2tt in the (x,(/>) plane, (a) two frequency quasiperiodic torus at e=0.03, 
(b) doubled torus attractor for e=0.0312, (c) wrinkled doubled attractor for e=0.0335, (d) strange 
nonchaotic attractor for e=0.0342, (e) strange nonchaotic attractor for e=0.0345, (f) wrinkled at- 
tractor for e=0. 03461 (g) doubled torus attractor for e=0.0347, (h) merged attractor for e=0.036, 
The other parameters are w„=0.25, A=0.5, a=0.2, u p =1.0, 0,^=6.7 and w e =0.991. 

FIG. 9. Largest Lyapunov exponent X ma x vs.e corresponding to Fig.8 
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FIG. 10. Projection of the two frequency quasiperiodic attractors of Eq.(5) for /=0.342: 
Poincare plot with <f> mod 2ir in the (x,<j>) plane, (a) two frequency quasiperiodic torus at e=0.03, 
(b) doubled torus attractor for e=0.0309, (c) wrinkled doubled attractor for e=0.033, (d) strange 
nonchaotic attractor for e=0.0337, (e) chaotic attractor for e=0.034, (f) strange nonchaotic attrac- 
tor for e=0.0345, (g) doubled torus attractor for e=0.0347, (h) merged attractor for e=0.036, The 
other parameters are cj„=0.25, A=0.5, a=0.2, oj p =1.0, 0,^=6.7 and oj e =0.991. 



FIG. 11. Largest Lyapunov exponent X max vs.e corresponding to Fig.10 



FIG. 12. Projection of the two frequency quasiperiodic attractors of Eq.(5) for /=0.345: 
Poincare plot with <j) mod 2ir in the (x,0) plane, (a) two frequency quasiperiodic torus at e=0.03, 
(b) doubled torus attractor for e=0.0305, (c) strange nonchaotic attractor for e=0.0334, (d) dou- 
bled strange nonchaotic attractor for e=0.0337, The other parameters are u^=0.25, A=0.5, a=0.2, 
w p =1.0, Ql=6.7 and w c =0.991. 
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